Complex RG flows for 2D nonlinear 0{N) sigma models 
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Motivated by recent attempts to find nontrivial infrared fixed points in 4-dimensional lattice gauge 
theories, we discuss the extension of the renormalization group (RG) transformations to complex 
coupling spaces for 0{N) models on L x L lattices, in the large- limit. We explain the Riemann 
sheet structure and singular points of the finite L mappings between the mass gap and the 't Hooft 
coupling. We argue that the Fisher's zeros appear on "strings" ending approximately near these 
singular points. We show that for the spherical model at finite A'' and L, the density of states is 
stripwise polynomial in the complex energy plane. We compare finite volume complex fiows obtained 
from the rescaling of the ultraviolet cutoff in the gap equation and from the two-lattice matching. 
In both cases, the fiows are channelled through the singular points and end at the strong coupling 
fixed points, however strong scheme dependence appear when the Compton wavelength of the mass 
gap is larger than the linear size of the system. We argue that the Fisher's zeros control the global 
properties of the complex fiows. We briefiy discuss the implications for perturbation theory, proofs of 
confinement and searches for nontrivial infrared fixed points in models beyond the standard model. 
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I. INTRODUCTION 

Recently, there has been a renewed interest [Il-Q the 
possibility [13) of finding nontrivial infrared (IR) fixed 
points in asymptotically free gauge theories with enough 
matter fields. A particularly interesting situation is when 
in addition of the nontrivial IR fixed point one ultraviolet 
(UV) fixed points also appears at larger coupling. It 
has been argued [ll|, that in this type situation, a 
parameter can sometimes be varied in such a way that 
these two fixed points coalesce and then disappear in the 
complex plane. 

This observation has motivated us [T3| to study com- 
plex extensions of renormalization group (RG) fiows in 
the complex coupling plane. The main result is that the 
Fisher's zeros - the zeros of the partition function in the 
complex coupling plane - act as a "gate" for the RG flows 
ending at the strongly coupled fixed point. This can be 
seen as a complex extension of the general picture pro- 
posed by Tomboulis jT4| to prove confinement: the gate 
stays open as the volume increases and RG flows starting 
in a complex neighborhood the UV fixed point (where we 
have asymptotic freedom) may reach the IR fixed point 
where conflnement and the existence of a mass gap are 
clearly present. 

More generally, constructing RG fiows in complex 
spaces could improve our understanding of the conver- 
gence of expansions (such as weak coupling and strong 
coupling expansions) that are used in the neighborhood 
of fixed points. Even though, complexification is often 
used for hydrodynamical flows, we are only aware of 
two previous studies of complex RG fiows: one for ex- 
actly solvable lattice models and one discussing the 
possibility of chaotic behavior in the decimation of one- 
dimensional Ising models with complex coupling [l6j . 



In the following we discuss two complex extensions of 
RG fiows for 0{N) models on Lx L lattices, in the large- 
N limit. The models are introduced in Sec. |lll We pro- 
vide a closed form expression for the partition function in 
the approximation where the non-zero modes of the La- 
grange multiplier are neglected. This is justified in the 
large- limit where we have equivalence with the spher- 
ical model. In Sec. IIII[ we study the map between the 
mass gap and the 't Hooft coupling A* — 1/6. We 
show that the map requires a Riemann surface with q + l 
sheets and 2q cuts in the A* plane, where q is an integer of 
order L^. By connecting the sheets in a specific way, we 
construct one circle at infinity in the A* plane (or around 
in the b plane) that maps into the circle at infinity in 
the plane and q others that maps in small regions 
near real interval [—8, 0]. 

In Sec. HYI we use the closed form of the partition 
function of Sec. |lI]to calculate the Fisher's zero at finite 
L and N. We show empirically that these zeros appear 
on "strings" coming from infinity in the b plane and end- 
ing near the singular points of the map discussed in Sec. 
mil This is consistent with the infinite volume picture 
provided in Ref. [l^]- The density of zeros on these 
strings scales like and N. The results of this section 
can be compared to what is found for other models, for 
instance in Refs. [T8l - [2l| where one and two-dimensional 
structures have been observed. 

In Sec. |Vl we show that the density of states is 
piecewise polynomial on q horizontal strips in the en- 
ergy plane. We discuss the conjecture that connects the 
Fisher's zeros with the zeros of the second derivative of 
the logarithm of the density of states. In Sec. I VII we ex- 
tend two RG methods to the complex b plane. The first 
one is based on a simple rescaling of the cutoff in the gap 
equation. The second one is a procedure called the two- 
lattice matching 0, [2^ . All the numerical calculations 
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are done with L and N even. 

Before embarking into the technical discussion, it is im- 
portant to keep in mind our various motivations. From 
a practical point of view, it is easier to calculate Fisher's 
zeros than to construct RG flows and establishing a 
clear connection should provide more robust ways to de- 
cide about the existence of nontrivial IR fixed points. 
The more general question of understanding the analytic 
properties of the map between the coupling and the mass 
gap in the complex plane is also important and includes 
the correction to asymptotic scaling. This question can 
be studied explicitly for the models considered here and 
the results illustrate the intricate pattern that can be 
produced by combining lattice artifacts and finite size ef- 
fects. Possible implications for lattice gauge theory will 
be discussed in the conclusions. 



For N even, the exponential of -V{N/2)2, (M^) is a 



product of poles located at the real negative values 



Mi 



-2^(l-cos(fc,)) 



(4) 



i=l 



The integer j indexes the various values taken while k 
runs over its V possible values k = ^n. The number of 
distinct poles will be discussed in section Hill where it is 
denoted q + 1. The number of times a given value of M'j 
occurs will be denoted rij and we have 'Y^- nj = V. 

By calculating the residues, we get a general expression 
of the form 



1 VNbM^ 



(5) 



II. THE MODEL 

The partition function for the 0{N) nonlinear sigma 
model on a square, or more generally hypercubic lattice 
with volume V = reads: 



Z = C l[d^<j,^Si$^.$^~l)e- 



(l/9o)E=c,e(l-'#'-'Ax + e) 



(1) 

We denote the inverse 't Hooft coupling as 6 = 1/A* = 

l/g^N. The constraint ^x-^x = 1 can be implemented 
by introducing a Lagrange multiplier at every lattice 
site. After this is done, the action becomes quadratic in (jj 
and the Gaussian integration can be performed. It can be 
shown that in the large- TV limit only the zero mode of the 
Lagrange multiplier, denoted hereafter, survives [g^]. 
With this simplification, the partition function becomes: 



where aij are coefficients depending on the order i and 
the pole Mj. For a given j, the pole is of order njN/2 
and the index i runs between VN/2 — 1 and {V — nj)N/2. 
An explicit expression for 13 = 2, L = A, N = 2 is given 
in Eq. (A.l) in the Appendix. Despite the apparent 
singularities at 6 = 0, Z (6) is an entire analytical function 
and has a regular expansion at b = 0. For instance for 
D = 2,L = 4:,N = 2 



Z{b) = 1 - 646 + 
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0(6^) (6) 



If in addition L is even, then for every Mj, there is an 
associated Mj, = —8 — Mj obtained by changing all the 
kj into TT — kj and one can see that a, 



guarantees that 
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m = 
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with 



^ (^^') = 7 E ln[2 E (1 - cos ih)) + M^] 



(3) 



The contour of integration, denoted C, encircles the real 
interval [-8,0], also called "the cut". When Reb > 0, C 
can be deformed into a vertical line with an arbitrary pos- 
itive real part and a semi-circle at infinity going counter- 
clockwise from 7r/2 to 37r/2 which gives no contribution 
in the limit of infinite radius. Similarly, when Reb < 0, 
C can be deformed into a vertical line with a negative 
real part smaller than -8 and a semi-circle at infinity go- 
ing clockwise from 7r/2 to —tt/2. The prefactor has been 
adjusted in order to have Z{0) = 1 and the finite volume 
momenta take the values k 
integers modulo L. 



with n a vector of 



as explained in [T?']. 

It should be noted that the number of independent 
grows like x N. This proliferation of terms makes 
calculations performed in the next sections slow when 
L OT N becomes too large. For illustrative purpose, we 
will often use L — 4 and N = 2. This allows us to 
give explicit formulas of decent size as in the Appendix. 
However, it should be kept in mind that Eq. (2) is only 
a good approximation of the original partition function 
(Eq. (1)) for large TV. 



III. THE GAP EQUATION, SINGULAR 
POINTS, AND CUTS 

In the large- limit, it is possible to calculate the parti- 
tion function in the saddle point approximation. Varying 
A'P, we obtain the gap equation: 



d£.{M^)/dM^ = <B(M2) 



(8) 
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with 



^iM')^{l/V)J2 



2(E^i(l-cos(fc,)) + M2 



(9) 



9+1 



Af 2 - Mf 



(10) 



where nj is the number of times the pole Mj appears 
in the sum over the k and q + 1 the number of distinct 
poles. The explicit form of !B(M2) in the case D = 2 
and L = 4 is given in Eq. (A. 2) in the appendix. From 
now on Z? = 2 is assumed in all the examples. The 
precise value of q depends on accidental degeneracies but 
generally increases like . From Fig. [1] and Table HH 
we see that most of the values of 2q follow the relation 
(L/2 + l)2 or (i/2 + l)2 — 1 with exceptions every three or 
six data points. In general, after reducing to a common 
denominator, we obtain a rational form : 

*B(Af^) Q{M^)/P{M^) , (11) 

where Q and P are polynomials of degrees q and q + 1 
respectively. 




FIG. 1. Relation between L and 2q for D = 2. 

We now discuss the poles, zeros and singular points of 
the mapping between h and given by the gap equa- 
tion ([8]). From Eq. (g]), the q + 1 poles of *B are real 
and between -8 and 0. -8 and are always poles and on 
the real interval between them (that we call "the cut" 
hereafter), ©(M^) is zero once between each pair of suc- 
cessive poles. In addition of these q zeros, 05(1/2) is also 
zero when becomes infinite. This also makes q-\-l ze- 
ros. In general, h = S(A/2) takes all the complex values 
q + 1 times when is varied over the whole complex 
plane. Thus, the inverse map between the mass gap 
and h requires a Riemann surface with q + 1 sheets in 
the h plane. To decide where to put the cuts and how to 
join different sheets, we need to study the singular points 
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FIG. 2. Zeros, poles and singular points of b{M^) in the 
plane for 4x4 and 8x8 lattices. 



where db/dM"^ = 0. This occurs when P'Q~PQ', a poly- 
nomial of degree 2q, vanishes. The 2q roots of P'Q — PQ' 
appear in complex conjugate pairs in plane. This is 
illustrated in Fig. [2] for L = 4 and 8. We notice that 
as L increases the region where the singular points ap- 
pear shrinks along the cut. A log-log plot of the largest 
imaginary part of the singular points versus L is rather 
irregular but suggests that the height of the region where 
the singular points appear is of the order 1/ L. In the in- 
finite volume limit, the singular points become dense and 
cover the [—8, 0] cut. 

The image of a singular point M^^^^ in the h plane is 
*B(Af2.^^ )^ At infinite volume, *8(M2) becomes an inte- 
gral with four logarithmic singularities [l7| . The image 
of two lines of points located very close above and below 
the [—8, 0] cut, span four curves forming a cross shaped 
figure that can be seen in Fig. [31 For comparison, the 288 
singular points for a 32 x 32 lattice are also displayed. 
We find that the real part of the closest singular points 
(CSP) move to infinity while the imaginary part stay at | 
as the volume increases. Near a singular point M'^^^g , we 
have («8(Af2^_^ + 5z) - S(M2„g )) oc {5zY and we need 
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two sheets to invert the function in the neighborhood of 
03 (M^ ). 
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FIG. 3. The blending small crosses (x, blue online) are the 
b images of two lines of points located very close above and 
below the [—8, 0] cut in infinite volume; the crosses (+) are 
the images of the singular points for L = 32. The images of 
the 4 closest singular points appear as boxes. 

In order to construct the 9 + 1 sheets, we start with 
the region of the plane where + 4| >> 8. In this 
region we have b ~ We call this sheet the "main" 

sheet because it contains the usual strong coupling re- 
gion where b is small, real and positive corresponding 
to a m'^ap large, real and positive. As we now consider 
smaller values of |M^-|-4|, and correspondingly larger val- 
ues of b on the main sheet, we start running into singular 
points and need to decide on the location of the cuts. A 
simple choice is to take the cuts on vertical lines in the b 
plane going from the images of singular points with posi- 
tive imaginary part to increasing values of the imaginary 
part and from the images of singular points with nega- 
tive imaginary part to decreasing values of the imaginary 
part. The cuts are shown for L = 4 in Fig. |4l We can 
now construct the inverse image of the two branches of 
a cut on the main sheet. They end up on two real nega- 
tive values of where b becomes infinite. The complex 
conjugate of the inverse image of these two branches cor- 
responds to the complex conjugated cut in the b plane. 
Joining the two together, we obtain an oval shaped re- 
gion in the plane located symmetrically across the 
cut. If we vary inside each of the oval shapes, b runs 
over the whole complex plane forming the other q sheets. 
This construction is illustrated for L = 4 where g = 4 in 
Fig. m If a curve in the plane enters an oval shaped 
region say on the left of the critical point and exits on the 
right of this critical point, then its image in the b plane 
will wrap around the image of the singular point. In Sec- 
tion VI, we will show that the cuts in the b plane and 
the boundaries of the oval shaped regions are important 
to understand the RG flows. 
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FIG. 4. Singular points and cuts in the 6 plane for L = 4 
(top) and their inverse images in the plane. 



IV. FISHER'S ZEROS 

In this section, we discuss the Fisher's zeros of the 
partition function at finite N and L and the way their 
density scales with these quantities. Later we will show 
that these zeros play an important role in controlling the 
RG flows. The coefficients Uij in the large- expression 
Eq. ([5]) can be calculated exactly at finite A^ and L using 
the residues theorem. We can then search for the zeros 
of the partition function by using Newton's method and 
check that the number of zeros found inside a given region 
of the b plane encircled by a closed curve C is consistent 
with 



^ db{dZ/db)/Z ^ i2TTj2Qnq{C) , 



(12) 



where nq{C) is the number of zeros of order q inside C. 
The results are shown in Fig. [5] for L = 6. We see 
that the zeros form linear structures ("strings") ending 
at locations close to the (A^- independent) singular points. 
Similar pictures are found for other not too large values 
of L and A^ where similar calculations are feasible. In 
all the examples considered, we also found that the zeros 
closest to the real axis always have an imaginary part 
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larger than 1/8 in absolute value, in other word they 
never get closer to the real axis than the CSP. 

The density of zeros increases with N and L. We calcu- 
lated the density of zeros in the b plane (number of zeros 
in a given area of the b plane taken as large as possible) . 
The results are shown in Fig. HI The fits of these log-log 
plots show that at fixed L = 2, the density grows like 
^1.000 j^^j TV = 2, the density grows like £2-027^ 

This data is consistent with the idea that the density of 
zeros increases like the number of fields {NV). 




FIG. 5. Zeros of partition function for L = 6, N — 2 (boxes), 
and images of the singular points of b{M'^) (crosses). The 
images of the solutions /" — discussed in section |V] are 
given with the third symbol 
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FIG. 6. Number of zeros in a fixed region of the 6-plane for 
N — 2, L variable (top) and L — 2, N variable (bottom) 



DENSITY OF STATE 



Another way to obtain information about the location 
of the Fisher's zeros is to calculate the density of states 
n{E) in the complex energy plane. First, we consider the 
case where the energy E is real. The density of states 
n{E) is the inverse Laplace transform of the partition 
function: 



n{E) 



N 



K-^ioo 



dbe'^'^'Zib) 



(13) 



The contour of integration is a vertical line in the complex 
b plane with a constant positive real part K otherwise 
arbitrary. For L even, the relation between Z(b) and 
Z{—b) given in Eq. (O implies that 



n{W -E)^ n{E) . 



(14) 



We can now use the form of the partition function given 
in Eq. ^ to obtain an explicit form. The only poles of 
the integrand are at 6 = 0. If E+ ^Mf > 0, we can close 



the contour by adding a semi-circle at infinity to the left 
and calculate the residue of the pole of order j at & = 0. 
If -I- < we can close the contour by adding a 

semi-circle at infinity to the right and the closed contour 
includes no poles. Since all the M| are real and negative, 
it is clear that n{E) ^ for E <0. For L even, Eq. 
implies that that n{E) = for i? > W. The final result 
for E real and L and N even is that the density of states 
is piecewise polynomial for < < AV and zero outside 
this interval. An explicit form will be given in Eq. (|15|). 



We now generalize this construction to the case where 
E is complex. If we consider e''-'^^ on a circle at infinity 
in the b plane, the expression blows up on one-half of the 
circle and decays on the other half. If we insist on being 
able to define the density of states by integrating Eq. ([5]) 
term by term then the only way to extend the definition 
is to rotate the line integral in such a way that bE is 
purely imaginary at both ends. The argument about the 
closing of the contour goes as before and we enclose the 
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poles at if Re{E + ^Mf) > 0. The final result is: 



L=2, N=4 



V 



(^-l)! 



e{Re{E+-M^)) 



(15) 

FigEl shows the function for L = 4 and N — 2 following 
an explicit formula given in Eq. (A. 3) in the appendix. 




FIG. 7. Density of state function for L = 4 and N = 2 

In summary, we have constructed a complex exten- 
sion of the density of states that is stripwise polyno- 
mial in the complex E plane. As the polynomials as- 
sociated with n{E) in two contiguous strips are dif- 
ferent, it is unavoidable that some of their derivatives 
will be different at the boundary. However the discrep- 
ancies only appear at some order that increases with 
VN/2. This can be seen from the discussion in Sec- 
tion ini where we discuss the range of powers of 1/6 
appearing in the partition function. As we cross the 
boundary of a strip Re{E) = V\M'j\/2, we add terms 
of the form [E -\- ^M'jy~^ which vanish at the bound- 
ary if i > 1. This term generates a discontinuity in 
the i — 1-th derivative. The lowest value of i occurring 
is iy — nj)N/2, where rij is the number of times the 
pole Mj appears. Consequently, the lowest derivative at 
which a discontinuity occurs is (V — MaXj{nj))N/2 — 1. 
For instance, for L — 4, Maxj{nj) — 6, and the low- 
est order is 5A^ — 1. It can indeed be proven that for L 
even, Maxj{nj) = 2{L — 1) and that it corresponds to 
M| = -4 for which "mirror" momenta can be paired in 
maximal number. 

Interestingly, we can use the large- iV limit to obtain 
finite volume thermodynamics from the density of states. 
At large N, by saddle point approximation of the Laplace 
transform of Eq. ([T3|) . we obtain that 



biE) EE f = 



n'{E) 



(16) 



Nn{E) ' 

which is similar to the standard thermodynamical rela- 
tion /3 = dS/dE. 
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FIG. 8. Singular points, zeros of partition function, and 
Re(f") = in £ plane for L = 2,N = 4. The line (red 
online), corresponds to Re{f") = with Re{f") > above 
it. 



Eq. (fTB|) provides a mapping of vertical strips of the 
E plane into the &-plane. The solutions of f"{E) — 
give the singular points of this mapping. These singular 
points can be found strip by strip by finding the solutions 
of the polynomial equation: 



n{E)n"{E) = in'{E)y 



(17) 



that belongs to that vertical strip (assuming that n{E) ^ 
0). We can then compare the 6-images of these singular 
points using Eq. strip by strip. Fig. [S] shows the 

images of these singular points for L = 6, = 2 together 
with the images of the singular points of 6(Af^) and the 
zeros of partition function in h plane. We see that the sin- 
gular points of h(E) seem to cluster is the same region as 
the singular points of 6(Af^) but with a slightly broader 
range. We conjecture that in the infinite volume limit, 
both types of singular points will accumulate on the out- 
side boundary of the crossed shaped curve displayed in 
Fig. [3] (image of the cut in the infinite volume limit). 

In the saddle point approximation, it can be argued 
|25j , that Fisher's zeros can only appear as images of 
regions in the E plane where Re{f") > 0. In Fig. |8l we 
see that the inverse image of the Fisher's zeros appears 
in the region Re{f") > even at small values of N and 
L. 



VI. COMPLEX RG FLOWS FOR 0{N) 

In this section, we construct complex RG flows using 
two different methods. The first one is based on a rescal- 
ing of the UV cutoff in 6(M^) given by the saddle point 
relation Eq. ([5]), the second is a complex extension of the 
two-lattice matching proposed in Ref. [l^, 
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A. Rescaling of 

Eq. (|5]) can be interpreted as a relation between the 
bare coupling and the UV cutoff keeping the renormalized 
mass, or mass gap, fixed. As the UV cutoff is lowered, the 
coupling increases and ultimately, b flows to zero. More 
speciflcally, Af^ = m|,/A^ and we follow the change in b 
under the change A — )■ A/s. In the literature the rescal- 
ing factor s is often denoted 6, but we are already using 
this symbol for the inverse 't Hooft coupling. Under this 
change, — >■ s^M^ and we can follow the trajectory in 
the b plane. In the following, we will iterate the transfor- 
mation with s = 2 in order to allow a comparison with 
the other method for which s = 2 is the simplest possibil- 
ity. The initial values of were taken on a small circle 
around the origin in the plane and then multiplied 
repeatedly by 4. We can visualize these trajectories as 
"rays" coming out of the origin in the plane. 

At infinite volume, as long as the trajectories in the 
do not cross the cut, the corresponding trajectories 
in the b plane will stay inside the cross shaped image of 
the cut shown in Fig. |31 Sample trajectories are shown 
on Fig. ini a) and illustrate this idea. The fiatness of 
the fiow at larger values of Reb can be understood from 
the approximate logarithmic scaling of b and the fact the 
rescaling factor is real and does not affect the phase of 
the rays. 

At finite volume, the cut acquires a thickness and a 
structure described in Sec. IIIIl If the ray crosses the 
oval shaped regions on two sides of the singular point 
(see Fig. 3]), then the corresponding trajectory in the b 
plane will wrap around the image of the singular point 
in the b plane using another Riemann sheet. At very 
small AP, we have b{AP) ~ 1/{L'^A'P) (instead of a 
logarithmic dependence) and we have inverted rays at 
infinity in the b plane. By small, we mean that the other 
poles at Alj can be neglected which occurs if AP << 
1 /i^, in other words, if the Compton wavelength is larger 
than the volume of the system. Sample trajectories are 
shown in Fig. [HI b) for L = 6. The two singular points 
which have the closest distance to the real axis are called 
the closest singular points (CSPs). 

B. Two- lattice matching 
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FIG. 9. a) (up) RG flows by rescaling and image of the cut 
at infinite volume; b) (bottom) flows by rescafing, singular 
points, CSPs, zeros of partition functions (A*' = 2), and /" = 
(A'^ = 2) in 6 plane for 6x6 lattice system, c) RG flows for 
the 2-lattice matching between 8x8 and 4x4 lattices. Circles 
and triangles are the singular points for L = 4 and L = 8. 



The 2-lattice matching [2^ l23| is another method that 
can be used to obtain complex RG flows. We will be 
matching observables i?(&, L) from systems with different 
lattice spacing but equal physical sizes. Under a RG 
transformation, the lattice spacing increases (a — sa) 
and the number of sites decreases (L — >■ i/s), but the 
physical length La stays constant. In the following, we 
compare a 2" x 2" lattice with 't Hooft coupling b with 
a 2"'i X 2"~i lattice with 't Hooft coupling 6'. The 
2" X 2" model is blocked n — 1 times while the coarser 
lattice model is blocked n — 2 times. For an arbitrary 
lattice with L = 2"?, we deflne the ratio of correlations of 



block observables (sum of all the spins inside a. L/2x L/2 
block B or its nearest neighbor block NB). 



R{b, L) 



eNB '^y)) 



b.L 



(18) 



Note that R{b, L) is independent of the rescaling of the 
flelds that needs to be performed in order to get an ex- 
plicit form for the RG transformation. Note also that 
in order to deflne the numerator and denominator sepa- 
rately, we need to divide by the partition function, but 
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these normalization factors cancel in the ratio. In the 
large- limit, the correlations are the same as for a Gaus- 
sian model with a mass defined as a function of b by 
Eq. ([5]) . By using the binary decomposition of the inte- 
gers between 1 and 2"^^, in the Fourier decomposition 
of the two-point function, we obtain 

( (E ^-)( E ^vn = E^°«(2"-ifci)jo(k)0(k) , 



x£B yeNB 



6,2" 



(19) 



xGB yes 



b,2" 



with 



Sj{k) - n + cos(2'A:i))(l + cos(2'fc2)) , (20) 
1=1 

and the lattice propagator 

e5(k) = l/[2(2 - cos(fci) - cos(fc2)) + M^] . (21) 

After performing the sums over the momenta, R{b^ 2") 
reduces to a ratio of polynomials in NP understood as 
function of b. Examples are given in Eqs. (A. 4-7) in 
the appendix. As explained in Sec. Illli there are q + 1 
values of corresponding to one values of b. In the fol- 
lowing, we always select the value corresponding to the 
"main sheet". This includes the circle at infinity where 
the conventional strong coupling behavior applies. In 
other words, we exclude the q other values of cor- 
responding to the small oval shaped region surrounding 
the cut. 

A complex RG map can be constructed as follows. 
Given an initial complex value of 6, we determine M^(6) 
corresponding to the main sheet. Using Eqs. (|18m9p . 
this result in a unique numerical value for 2"). We 
then match this number with R{b' ,2'^~^) expressed as a 
ratio of polynomials in M' . This results in a certain 
number of solutions for M'^ . We only keep the ones 
corresponding to the main sheet. Each of these selected 
solutions determines a unique value of b' . If more than 
one remain, we only keep the one closest to b. In order 
to quantify the level of ambiguity associated with this 
choice, we define the ambiguity A as 



A 



16-6' 



|6-6'|„ 



ext — to—min 



(22) 



where |6 — b'\min is the distance between 6 and the clos- 
est solutions and |6 — b'\next-to-min the distance to the 
next closest solution. If there is only one solution ^ = 
(|6 - b'\next-to-m.in = oo). There is also the logical pos- 
sibility that none of the M'^ solutions correspond to the 
main sheet, but we never encountered this case in practi- 
cal calculations. In Fig. [3};), light regions correspond 
to small values of A while dark-colored regions stand 
for values of A close to 1 (maximal ambiguity). When 



|6-6'|m„ ~ \b-b'\next-to-mm it is uot clcar that we can 
define a RG flow, and not surprisingly chaotic behavior 
is often observed in these circumstances. Unambiguous 
flows tend to stay in light-colored regions while ambigu- 
ous flows go into the dark regions and jump to other fixed 
points. By looking at all the singular points, we find that 
the closest singular points (CSP) are close to the edges of 
dark-colored regions, which gives us the idea that the un- 
ambiguous RG flows are bounded by the singular points 
and the features of all the RG flows are controlled by the 
CSP. Fig. ini shows the RG flows starting from the first 
fixed point on the real positive axis(6o ~ 0.64) with the 
singular points of 4 x 4 and 8x8 lattice systems. 

In order to understand the fixed points on the real axis 
more systematically, we study the RG transformation in 
the small M^, large 6 limit. For very small M^, the pole 
at zero dominates in Eq. ([8]) and we have 6 ~ 1/(M^L^) 
. The other poles have small contributions provided that 
« 1/L2^ in other words, when the Compton wave- 
length is larger than the linear size of the system. In the 
same limit, the matching condition becomes 



M' ~ 4(1 - A/L^)M^ 



(23) 



with A a constant that we will determine later. We now 
only consider the real solutions of 6 and calculate the 
change of the coupling A6 = 6 — 6'. Numerical results are 
shown in Fig. 1101 where a comparison with the rescaling 
method is made. In the infinite volume limit, the rescal- 
ing from to AM"^ leads to the relation: A6 = 77- In 2 
and for large 6, there is no nontrivial fixed point. Putting 
everything together, we get the approximate finite vol- 
ume formula: 



A6(6)^-(A/i>+i^ 



(24) 



This implies the that we have an approximate fixed point 



b* 



In 2 

A2^' 



(25) 



In Table m we compare the numerical fixed points 6*„„ 
with the approximate b*pp with A — 30.5, and find 
the approximate model provides reasonable estimates for 
large volume. For comparison, A6 can also be calculated 
with the rescaling method. In this case M' = 4M2 and if 
A/^ is finite and nonzero, the denominator of all the terms 
is 6' are strictly larger than for 6, consequently 6' < 6 and 
there is no nontrivial fixed point. In summary, we see 
that discrepancies between the two RG methods occur in 
the limit where the Compton wavelength is larger than 
the size of the system. Otherwise, the RG flows going 
through the CSP are very similar for the two methods. 



C. Comparison of the two methods 

In this subsection we argue that in the infinite volume 
limit, the two RG fiows discussed in the two previous 
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L 




^app. 


4 


0.320 


0.058 


8 


0.648 


0.23 


16 


1.47 


0.93 


32 


4.36 


3.70 


64 


15.5 


14.82 


128 


59.9 


59.3 


256 


237 


237 



TABLE I. Numerical solutions of fixed points and approxi- 
mate solutions from Ab = 



subsections should coincide. By construction, the infinite 
volume of the RG flows illustrated in Fig. [H]b) turn into 
those of Fig. |9]a). It is nevertheless interesting to figure 
out in detail how it occurs. As we take initial conditions 
corresponding to = ee*^, we see that as 9 becomes 
slightly larger than tt/2, the linear RG trajectories in the 
plane cross the riglitest oval shaped region discussed 
in Sec. Illll on the both sides of the upper singular point. 
As a consequence, the corresponding RG flows in the 6- 
plane wrap around the CSP which is the image of the first 
singular point, before going to 0. In the infinite volume 
limit, the CSP moves to infinity and we recover Fig. ^a). 
In some sense, the finite volume provides a regularization 
and gives a mathematical meaning to what happens for 
initial conditions corresponding to the cut. 

The two-lattice matching relies on the fact that if we 
apply the RG transformation enough times, the RG flows 
are projected on the unstable manifold which is one- 
dimensional for the models considered here. In this limit, 
the flows from the two lattices can be compared unam- 
biguously. From the way we set up the calculation, it is 
clear that we reach an infinite number of RG transfor- 
mation in the infinite volume limit (we block-spin ln2 (L) 
times). In this limit, the finite volume fixed point goes to 
infinity like and it is plausible that the unambiguous 
flows of Fig. [9]c) become similar to those of Fig. [9] a). 



VII. CONCLUSIONS 

In summary, we have extended two types of RG trans- 
formations to complex coupling spaces for 0{N) models 
on LxL lattices in the large- limit. The three graphs of 
Fig. [5] illustrate possible outcomes in models where cal- 
culations are more difhcult such as asymptotically free 
lattice gauge theories. Our general expectation is that at 
infinite volume, almost horizontal flows come from the 
large j3 region where the logarithmic scaling character- 
istic of asymptotic freedom holds. The flows are then 
funneled through a continuous boundary whose shape is 
a lattice artifact. For instance, if we had used a contin- 
uum model with a sharp cutoff, the boundary would be 
flat instead of the shape seen in Fig. [S] a) . The existence 
of this boundary reflects the fact that we try to extend 
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. 11 
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4 to 21 

16 to 8 
64 to 32 
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log (2) /2/Pi 
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FIG. 10. A6 versus b from rescaling (top) and 2-lattice match- 
ing (bottom). 



the mapping between and b, we encounter a cut in 
the plane and the boundary is the image in the b 
plane of a close curve tightly enclosing the cut. 

At finite volume, the continuous boundary is replaced 
by a loosely defined region where complicated or ambigu- 
ous trajectories are observed (see also Fig. 1 in Ref. 13]). 
Empirically, this region seems to coincide with the region 
where the images of the singular points of the mappings 
b{M'^) or b{E), and the ends of strings of Fisher's zeros 
appear. As the volume increases, the number of zeros in 
a fixed area of the 6-plane increases like the volume and 
N and we believe that the Fisher's zeros become dense 
outside of the boundary mentioned above. Fig. 3 in Ref. 
[l3t suggests that it will also be the case in lattice gauge 
theory. Our numerical study supports the idea that by 
monitoring the lowest zeros of asymptotically free theo- 
ries when the volume increases, we can determine if the 
RG flows reach the region where a mass gap is present 
or if instead a nontrivial IR fixed point is encountered. 
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We are planning to investigate the scaling of the Fisher's 
zeros in SU{3) Nj flavors and also to investigate the ques- 
tion of the corrections to asymptotic scaling [26| in the 
complex plane. 

The complex flows described here are in some sense 
the simplest possible ones and we expect similar complex 
flows for SU{2) lattice gauge theory in 4 dimensions. In 
contrast, for U{1) lattice gauge theory in 4 dimensions 
it ap pea rs that the Fisher's zeros pinch the real axis as 
[23] with 2 < X < 4. A more precise estimate for x 
will be discussed in a forthcoming preprint [28]. In this 
case, we expect complex RG flows similar to those of the 
3 dimensional Ising or 0{N) models where two phases 
are present. 

Complex flows for models with a conventional second- 
order phase transition have been constructed for the hi- 
erarchical model [13]. More recently, it has been found 
that the qualitative behavior of the complex flows can 
be modifled by lowering the adjustable parameter (usu- 
ally denoted c) below the critical value (c=l) where a 
second order phase transition is possible. As c reaches 
1, the nontrivial fixed point moves to infinity. As c is 
further lowered, a pair of complex conjugated nontrivial 
complex fixed points appear [29|. In general, it seems 
likely that by considering actions with tunable parame- 
ters, it is possible to create interesting patterns for the 
complex RG flows. For hierarchical models with slightly 
modifled interactions, it is possible to create additional 
flxed points and study their effects on the discrete A/3 
function. In particular, we would like to understand if 
by constructing a continuous Calan-Symanzik /3 repro- 
ducing approximately the discrete one, it is possible to 
relate the nontrivial complex flxed points to the complex 
zeros of the continuous /3 function. This will also be dis- 
cussed in Ref. 

For the models considered here, a richer complex flow 
behavior could be obtained by adding new terms in the 
energy function and considering higher dimensional com- 
plex flows. One possibility that comes to mind is to add 
a term inspired by the Witten-Wess-Zumino term in the 
continuum. Multidimensional RG flows can exhibit in- 



tricate global properties sorted in Ref. [s^. Their com- 
plexiflcation appear to be a completely open fleld of in- 
vestigation. 

ACKNOWLEDGMENTS 

Part of this work started during the workshop "New 
applications of the renormalization group method in nu- 
clear, particle and condensed matter physics" held at the 
Institute for Nuclear Theory, University of Washington, 
Seattle (INT-10-45W). Additional work was done at the 
Aspen Center for Physics in May and June 2010 during 
the workshop "Critical Behavior of Lattice models" . We 
thank the participants of these two workshops for stim- 
ulating discussions. This research was supported in part 
by the Depar tment of Energy under Contra ct No. FG02- 



L 


2q 
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2q 


L 


2q 


L 


2q 


2 


4 


16 


80 


30 


220 


44 


528 


4 


8 


18 


100 


32 


288 


46 


576 


6 


16 


20 


120 


34 


324 


48 


580 


8 


24 


22 


144 


36 


328 


50 


676 


10 


36 


24 


148 


38 


400 


52 


728 


12 


40 


26 


196 


40 


440 


54 


784 


14 


64 


28 


224 


42 


472 


56 


840 



TABLE II. L and 2q 



91ER40664. 



Appendix: Numerical examples 

The numerical values of 2q for L up to 42 are given 
in Table HIl We also provide the explicit form for the 
partition function, *B(Af^) and the density of states for 
D = 2, L = 4:, N = 2. 



1 '=i76'i7'i<=~i28b 

- 5e^2^^ + 5 + 40e^^'' (2949126^ - 2027526^ + 534726 - 5255) 
-I- 40 (294912&3 -I- 2027526^ -I- 534726 -1- 5255) e^^''] 

^„ 353286^ 232652863 7376076806^ ^ , 

= 1 - 646h \ \-0(W) (A.l) 

17 51 969 ^ ^ ' 



*B4x4(m2) = ^(l/M^ + 4/(2 + M^) + 6/(4 + M^) 
16 V 

+4/(6 + M^) + 1/(8 + Af2)) . (A.2) 
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n(E) = [-644087808 (£; - 48)"e'(£' - 48) 

^ ' 7968993439842526298112^ ^ ' ^ ' 

- 73801728(£: - iS^^eiE - 48) - 2994432(£; - 48)"6i(i; - 48) 

- 42040(i; - 48)"6'(i; - 48) - 10882507603968(i; - 32)^9{E - 32) 

- 32848478208(i; - 32)"6'(£; - 32) - 12845952(i; - 32)"6i(i; - 32) 

- 644087808(i; - 16)"6i(i; - 16) + 73801728(^ - 16)^^e'(£; - 16) 

- 2994432(i; - 16)^^6'(£; - 16) + 42040(i; - 16)"6'(i; - 16) 

+ E^^0{E) -{E~ 64)^'^e{E - 64)] . (A.3) 
We also give the explicit form of the first R{b, 2') with understood as a function of b as explained in Sec. I VII 

2 -I- A,f 2 

^(^'^)- 2 + 4M2 + (M2)2 ' (^-^^ 



and 



32 + 82M^ + MjMY + ISjM^f + {M^ 
^ ' '~ 32 + 222A/2 + 314(M2)2 + 153(Af2)3 + 30(A/2)4 + 2(M2)5 ' ^ ' ' 
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